Glass-like low frequency ac response of ZrBi2 and Nb single crystals in the surface 

superconducting state 
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We report experimental studies of the low frequency electrodynamics of ZrBi2 and Nb single 
crystals. AC susceptibility at frequencies 3 - 1000 Hz have been measured under a dc magnetic 
field, Hq, applied parallel to the sample surface. In the surface superconducting state for several Hq 
the real part of the ac magnetic susceptibility exhibits a logarithmic frequency dependence as for 
spin-glass systems. Kramers-Kronig analysis of the experimental data, shows large losses at ultra 
low frequencies (< 3 Hz). The wave function slope at the surface was found. The linear response of 
the order parameter to the ac excitation was extracted from the experimental data. 

PACS numbers: 74.25.Nf; 74.60.Ec 



I. INTRODUCTION 



Nucleation of the superconducting phase in a thin sur- 
face sheath in a decreasing magnetic field parallel to 
the sample surface was predicted by Saint- James and dc 
Gennes some years ago [l| . They showed that nucleation 
of the superconducting phase occurs in a magnetic field 
Hq < Hc3 = 2.39 kHc, where He is the thermodynamic 
critical field and k is the Ginzburg-Landau (GL) param- 
eter. Experimental confirmations of thisprediction were 
done in numbers of publications H, H, ^ IE 111 - It was 
found that in the in surface superconducting state (SSS) 
at low frequencies, the superconducting compounds are 
lossy materials. However, the losses in the SSS, which 
exceed the losses of the normal state, shows a peak at 
external fields Hq, when Hc2 < Hq < Hc3- This peak 
exists even at frequencies of a few Hertz. In general, 
the response is nonlinear, and frequency dependent. The 
critical state model was found adequate for the descrip- 
tion of the experimental data |^. 

In the last few years the SSS has attracted renewed in- 
terest from various directions ^,Jj,9, 10, 11, 14, 15]. The 
stochastic resonance phenomena in the SSS for Nb single- 
crystal, were observed in the nonlinear low-frequency re- 
sponse to ac fields, 0. In Q it was assumed that 
at Ho — Hc3 the sample surface consists of many dis- 
connected superconducting clusters and subsequently the 
percolation transition takes place at H^^ — 0.81i?c3- The 
aramagnetic Meissner effect is also related to the SSS 
Voltage noise and surface current fluctuations in Nb 
in the SSS has been investigated |0. Surface super- 
conducting states were detected also in single crystals of 
MgB2 ^3 and ZrBi2 A new theoretical approach, 

based on a generalized form of the GL functional, was 
developed in [Tlj . Surface superconductivity in three di- 
mensions was studied theoretically in 12, 13J. It was 
found 0,01 that, in general, the wave form of the sur- 
face current in an ac magnetic field has a non-sinusoidal 



character. The latter one can be described by a simple 
phenomenological relaxation equation for transitions be- 
tween metastable states . The relaxation time in this 
equation depends on the deviation from equilibrium and 
increases with decreasing of the excitation frequency |l5j . 

In spite of the extensive study, the origin of low fre- 
quency losses in SSS under weak ac fields, is not clear yet. 
The critical state model implies that if the amplitude of 
the ac field is smaller than some critical value, the losses 
are absent. Indeed, the experimental results for Pb-2%In 
alloy confirm this prediction !5|. On the other hand, the 
observed response 8J for an excitation amplitude of 0.01 
Oe (that, is considerably smaller than that used in |^), 
shows losses in SSS of Nb at 10 Hz. Our measurement on 
Nb single crystal at 733 Hz also has shown that the out- 
of-phase part of the ac susceptibility, x i is independent 
of the amplitude, and is finite at low amplitude values. 
We consider these results as confirming that the losses in 
SSS are caused by essentially linear phenomena and thus 
the critical state model cannot be used for an adequate 
description of the ac response. 

In this paper we present detailed experimental study of 
the linear low-frequency response of ZrBi2 {Tc = 6.06 K, 
K « 0.75) and Nb (T^ = 9.2 K, k w 1.5) single crystals 
in the SSS, for frequencies 3 < w/27r < 1000 Hz. We 
show that already at 3 Hz the ac susceptibility does not 
equal dM/dHo, where M is the dc equilibrium magne- 
tization. For some dc magnetic fields x oc Inw, which 
resembles spin-glass systems ^El- The Kramers-Kronig 
analysis of our experimental data predicts that for sev- 
eral dc magnetic fields huge loss peak should exist at very 
low frequencies. It is believed that the observed response 
presents the average value over many clusters, each of 
which is governed by a second order differential equation, 
with individual relaxation parameters. This is unlike a 
spin-glass system where the first order differential equa- 
tion can be used jlTj . The observed ratio of Hcs / Hc2 for 
our crystals differs from the predicted value - 1.69. It 
is believed that this is due to the nonzero slope of the 
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wave function at the surface, b = — -j-^f-- The increased 
ratio of Hj^Hc2 shows that the sign of h is unexpectedly 
negative p^ . 

The order parameter, 5", in the SSS has the form: 

y) = f{x)exp{iKy), (1) 

where f{x) equals zero inside the bulk of the sample {x- 
axis is assumed normal to the sample surface, Hq parallel 
to z-axis) and K is some as yet undetermined constant. 
In general case for any dc magnetic field there is a band 
of possible values for K for which the solution of the 
GL equations can be found. These solutions describe the 
metastable surface states with nonzero total surface cur- 
rent and only the solution with zero current corresponds 
to the equilibrium state. In this state the dc magnetic 
moment equals zero for H > Hc2- On the other hand, in 
the presence of an ac field, the nonzero response in SSS 
shows that the ac moment is not zero, thus the system is 
not in equilibrium. In the linear approximation 

K{t) = Ko + Rc[G{uj, Ho)h{Lu)exp{~iLut)], (2) 

where the external magnetic field is H{t) = Hq + 
Rc[ft,(a;)exp(— iwt)]. The function G{lu,Hq) characterizes 
the response of the order parameter to an ac field and is 
determined in this paper. 

II. EXPERIMENTAL DETAILS 

The measurements were carried out on ZrBi2 and Nb 
single crystals, which were grown in the Institute for 
Problems of Materials Science NAS, Ukraine, and in 
the Institute of Solid State Physics RAS, Russia, respec- 
tively. The dimensions of the crystals are 10.3 x 3.2 x 1.2 
mm'^ for ZrBi2 and 10 x 3 x 1 mm'^ for Nb. The details of 
the sample preparation and their magnetic characteris- 
tics were published previously V, 'l9] . All the magnetiza- 
tion curves were measured using a SQUID magnetome- 
ter. Ac susceptibility, in-phasc, x i E^^id out-of-phase, 
X , components were measured using the pick-up coils 
method 5, 20J. Each sample was inserted into one of a 
balanced pair coils. The unbalanced signal as a function 
of the external parameters: temperature, dc magnetic 
field, frequency and amplitude of excitation, was mea- 
sured by a lock-in amplifier. The experiment was carried 
out as follows. The crystal was cooled down at zero mag- 
netic field (ZFC). Then the magnetic field was applied. 
The amplitude and the phase of the unbalanced signal 
were measured in a given magnetic field, including zero 
field, at all frequencies. The amplitude of excitation was 
0.01-^0.5 Oe. We have supposed that in zero dc magnetic 
field the ac crystal susceptibility is equal to the dc sus- 
ceptibility in the Meissner state with negligible losses. It 
permits us to find an absolute value of the in-phase and 
out-of-phase components of the ac magnetic susceptibil- 
ity for all applied fields and frequencies. A "home-made" 
measurement cell of the experimental setup was adapted 



to a commercial SQUID magnetometer. The block dia- 
gram of the experimental setup has been published else- 
where jTHj . 

For the Fourier component of the magnetization 
m{t) = Iie[m(Lu, Hq, h(Lu))exp(—iujt)] in the linear ap- 
proximation one can write 

m{uj,Ho,h{uj)) = XniHuj) +ATTJs{uj,Ho,h{uj))/c)+ 
AnXsJsiuJ^Ho, h[uj))/c, 

(3) 

where Xn is the susceptibility of the normal core 
of the sample, and Js{u!, Hq, hitu)) is the Fourier 
component of the surface supercurrent, Js(i) = 
He{ Js{uj, Hq, h{uj)))exjp{—iLL!t). This equation takes into 
account the magnetic moment of the normal core of the 
bulk and the magnetic moment of the surface super- 
current 0. With Eq. © we can find Js{uj, Hq, h{uj)) 
and the surface susceptibility defined as Xs(ci;,-ffo) = 
Js{w, Hq, h(uj))/h{u!). The measured susceptibility , 
x{u!,Hq) = m{uj,HQ,h{u!))/h{uj) and surface susceptibil- 
ity, Xs, are connected to each other as follows: 

Xs{u^,Hq) = [x{u;,HQ)-x{u;,Hr,)]/[l+47Tx{io,Hr,)], (4) 

where Hn > Hcs is the magnetic field at which the sam- 
ple is in the normal state. These quantities, ( Js and 
Xs), characterize the response of the surface current and 
eliminate the contribution of the normal core in the bulk 
and the small unbalanced signal of the empty coils 

III. EXPERIMENTAL RESULTS 

The field dependencies of the ac susceptibility x , X of 
Nb single crystal at different excitation amplitudes (/iq) 
are presented in Fig. ^ The data show clearly that at 
the low amplitudes x , X are almost independent on the 
Hq. Similar results were obtained for the ZrBi2 crystal. 
We consider that at Hq < 0.05 Oe the observed response 
in both Nb and ZrBi2 crystals has a linear origin. 

Figs. E^i and Fig. ^ show the x , X and the ZFC dc 
susceptibility, Xdc = M/Hq, as a function of the dc field 
{Hq) of Nb and ZrBi2 respectively. The insets in Figs.|2t' 
andlSb present the ZFC isothermal magnetization curve. 
Note, that for both crystals, at Hq > Hc2 the dc suscepti- 
bility in the SSS is zero, whereas a large and diamagnetic 
ac susceptibility signal is observed, see Fig. El The dif- 
ference between the dc and the in-phase ac susceptibility 
signals remain even at a;/27r — 3 Hz (the lowest frequency 
in our set up). The out-of-phase component, x , has a 
broad maximum for Hq > Hc2- 

In the normal state, for Hq > Hcs, the losses for Nb are 
greater than for ZrBi2 (see Fig. 2a for 1065 Hz), because 
the normal state conductivity of Nb is larger than that of 
ZrBi2. The out-of-phase susceptibility X is proportional 
to Lu as expected for the normal skin effect in the limit 
6 >> d, where 6 is the skin depth and d is the sample 
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FIG. 1: (Color online) Magnetic field dependencies of x and 
X for Nb at T = 8.5 K at different amplitudes of excitation, 
ho- 



thickness ^J. For both crystals, when Hq > Hc3, x" is 
field independent. However in SSS, (see Fig.EJ, the losses 
for Nb are lower than for ZrBi2. Both, the in- and out-of- 
phase components of the ac susceptibility, xi^jHo), do 
not show any clear peculiarities as the dc magnetic field 
passes through Hc3- Therefore, it is difficult to extract 
the Hc3 value from the experimental data. More sensitive 
measurement of Hcs can be done by using the Xs{'-^i Hq) 
curves (see below). 

Figs.|31and0]demonstrate the frequency dependence of 
X and X for both Nb and ZrBi2 crystals at dc magnetic 
field above Hc2- The character of these curves depends 
on the field. Thus x is a concave function of lo at fields 
close to Hc2 and is a convex function close to _ffc3- For 
some fields, x cx: Inw for both samples. Such a logarith- 
mical dependence is typical of spin glasses which are com- 
pletely different physical systems . Fig. |31 presents 
X as a function of the reduced magnetic field, Hq/Hc2 
for two frequencies and temperatures for Nb and ZrBi2 
crystals. The Hc3/Hc2 ratio is temperature dependent 
as previously reported for Nb (e.g. |22j). 

IV. THEORETICAL BACKGROUND 

For interpretation of the experimental data we used 
the numerical approach to the normalized stationary GL 
equations 

-(iV/K+ A)2^'2 + ^_ I ^ |2 ^ ^ 

> > n (5) 

-curlcurlA = A | * ^ +i/2K(**V* - W^-*) 

for an external magnetic field parallel to the sample 
surface. The order parameter, is normalized with re- 
spect to its value at zero magnetic field, the distances 
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FIG. 2: (Color online) (a) Magnetic field dependencies of x i 
x" and Xdc = M/Ho for Nb at T = 8.5 K. 
(b) Magnetic field dependencies of x', x" ^ind Xdc = M/Hq 
for ZrBi2 at T = 4.5 K. 



with respect to the London penetration length A, and 

the vector potential A with respect to ^/2Hc\, where 
He is the thermodynamic critical field. Assuming that 
the order parameter has the form of Eq. |^ with the 
yet as undetermined parameter K, we see that K is the 
integral constant of these equations, if the sample thick- 
ness considerably exceeds the coherence length and the 
superconductor is homogeneous. The relaxation time in 
the time-dependent version of GL theory is of the order 
of 10~^^ — 10""'^^ s and we can use the stationary ver- 
sion of the GL theory, Eq. for the ac experiments. 
Equations |(SJ with proper boundary conditions and the 
requirement that the order parameter differs from zero 
only near the surface, have the solution for a whole band 
of K values. But only one value of K corresponds to 
the total surface current equals zero and this K de- 
scribes the equilibrium state with the minimal free en- 
ergy. These nonlinear equilibrium surface solutions have 
been discussed in detail in ,23j . The ac response for su- 
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FIG. 3: (Color online) Frequency dependence of the ac sus- 
ceptibility for Nb at T = 8.5 K at dc magnetic field Hq > Hc2- 
(a) - x'i^) and (b) - x"('^)- 



FIG. 4: (Color online) Frequency dependence of the ac suscep- 
tibility for ZrBi2 at T = 8.5 K at dc magnetic field Ho > Hc2- 
(a) - x'(t^) and (b) - x"('^)- 



perconductors in fields Hc2 < Hq < Hcs differs from the 
normal one. This means that in an ac field the super- 
conductor is in a nonequilibrium state with finite sur- 
face current and nonequilibrium K. Generally speaking, 
the total surface current depends on both: the instant 
values of external magnetic field H{t) and on ivr(t), thus 
Js{t) = Js{H{t), K{t)). If a small ac magnetic field h{t) 
is superimposed upon a dc field TJq, the amplitude of the 
surface current in the linear approximation is 



^-G{LJ,Ho)h{Lj). 



(6) 



dKo 



where Kq is the equilibrium value of if in a dc magnetic 
field Hq, and G(w, Hq) describes the linear response of K 
to an ac field Eq. The partial derivatives in Eq. © 

(4 ^ ''fj^f"^ and 4 ^ ^^'^smr^) can be calcu- 
lated numerically and in Fig.|Slwe show the results for Nb 
(GL parameter k = 1.5 ,22,j) and ZrBia {k = 0.75 pi). 

Equation ^ describes the linear response of the sur- 
face current to an external ac magnetic field. If K does 
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FIG. 5: (Color online) The field dependence of x for Nb 
at T = 4.5 K, 8.5 K and for ZrBi2 at T = 4.5 K (where 
uj/2tt = 765 Hz). 
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FIG. 6: (Color online) The derivatives of surface current Ju 
Jjf versus reduced magnetic field for Nb and ZrBi2 (see text). 



not change during an ac cycle, we see from Fig. 1^1 that 
only a smooth decrease of the surface current without 
any losses should be observed, as the dc magnetic field 
increases. The critical current model assumes that for 
h{uj) < he, where he is some critical ac field, the surface 
current follows the external magnetic field without any 
delay, ImG'(w, Ha) = 0, and therefore surface losses are 
absent. 
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V. DISCUSSION 

It is clear that the experimental data show the exis- 
tence of SSS in our crystals. Usually i/c3 is defined by 
the onset of the surface screening, i. e. by the appear- 
ance of the deviation of the ac response from its normal 
value, as the dc field decreases from some large value. Be- 
cause in the SSS the ac response is frequency dependent, 
using of low frequencies could give the underestimated 
value of Hcz- As an example, in Fig. [7^ we show a set 
of data obtained for Nb sample at 8.5 K, at frequencies 
115, and 765 Hz. The imaginary part of surface cur- 
rent is more sensitive for Hcz determination (obviously, 
the experimental setup sensitivity is higher at higher fre- 
quencies). The data at a frequency of about 100 Hz yield 
iJc3 = 680 Oe, while at higher frequencies (« 700 Hz) we 
obtain Hc^ = 760 Oe, in the later case Hc3/Hc2 = 1-9. 
This value is considerably larger than the value predicted 
1.69 in Ref. 1] With decreasing the temperature, the dis- 
crepancy between the theoretical prediction [ij and the 
experimental values increases and the ratio at T = 4.5 K 
becomes Hc3/Hc2 = 2.34 (Fig.[7|3). The decrease of this 
ratio for temperatures in the vicinity of Tc, was found in 
several experiments [22, |^ and it was associated with 
the decrease of Tc near the surface [2^ . In the framework 
of the GL theory the Hcs / Hc2 ratio can be changed only 



FIG. 7: (Color online) The field dependence of the real, J^, 
and imaginary, , parts of surface current at difi'erent fre- 
quencies for Nb at T=8.5 K , Hc2 = 400 Oe (a); 4.5 K, 
Hc2 = 3200 Oe (b). 



if (i) either Tc in the surface layer differs from the bulk 
value, or (ii) the slope of the wave function at the surface 
6 = differs from zero [ll|2|. In Fig. we show 

the calculated Hc3/Hc2 {k = 1.5) for these two cases. 

The boundary condition for the first case is 5 = but 
the Tc values of a surface layer with thickness L, increases 
as: 

(T,(x) - Tc)/iTc - T) - dTcexp{-x/L) (7) 

(Fig. IHli). In the second one, (when dTc — 0), we assume 
that b^O (Fig.Et). The ratio of Hc3/Hc2 ^ 1-9 at 8.5 
K, this can be a result of either the enhanced Tc by 0.13 
K (dTc = 0.2) at the surface layer with the thickness of 
L/X ~ 0.71, or, by assuming b — — 0.15/A. Decreasing 
the temperature results in an increase of both: L/X and 
the Hc3/Hc2 ratio. At T = 4.5 K, dTc = 0.028 and from 
Fig. one finds that the value of the Hc3 / i?c2 ratio 
cannot exceed 1.75. Therefore, if the GL theory is ap- 
plicable, than, at r = 4.5 K, the growth of the Hc3/Hc2 
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FIG. 8: (Color online) The Hci/Hc2 ratio for two cases: 

(a) Tc in the surface layer (of thickness L) is higher than in 
the bulk value see Eq. 0); 

(b) the slope of the wave function, h, at the surface is different 
from zero (k = 1.5). 



FIG. 9: (Color online) Real and imaginary parts of surface 
susceptibilities versus the dc magnetic field for Nb at 20 Hz 
and T = 8.5 K near the point where the absolute value of 
both components are equal to each other [Ho/Hca ~ 1-3). 



ratio to 2.34 in Nb crystal is due to 6 7^ 0, and the abso- 
lute value of b increases with decreasing temperature. 

The behavior of the ZrBi2 crystal is similar. The 
Hci/ Hc2 ratio increases from 1.5 at 5.5 K to 1.75 at 4.5 K. 
Since the ratio at 5.5 K is smaller than 1.69, it is pos- 
sible that Tc of the surface layer is smaller than that of 
the bulk. Our calculation shows that this value (1.5) can 
be obtained, for example, if dTc — 0.16 and L/\ — 2. 
For T = 4.5 K the dimensionless dT^ is equals to 0.05 
and Hcd./Hc2 ratio does not exceed 1.65. We notice that 
in order to explain the Hc3/Hc2 ratio, one has to take 
into account the nonzero slope of the wave function at 
the surface. 

In general, with decreasing the dc field from its maxi- 
mal value, both the real and imaginary parts of the sur- 
face current appear simultaneously, but at the beginning 
the imaginary part increases faster than the real one. 
On the other hand, at Hc2 a complete screening takes 
place and the absolute value of Xs reaches its maximal 
value, l/47r, while Xs = 0- So at some dc magnetic field 
I Xs I will be equal to | Xs I ■ This point was identified in 
Ref. !^ as the percolation transition from noncoherent 
SSS to the coherent one in the Nb sample. For our Nb 
crystal, at T = 8.5 K and frequency 20 Hz it occurs at 
Ho/Hc3 ~ 0.68, a value which is slightly smaller than the 
value 0.81 reported in Q. However, our data presented 
in Fig. |51 does not permit us to consider this point as a 
phase transition, due to a smooth maximum \ Xs I at this 
field and the absence of any peculiarity in Xs- 

The response of K to an ac field is described by the 
function G{uj,Hq) which can be found from Eq. 
Figs. [TUk. [TUb show G{uj, Hq) as a function of a reduced 
magnetic field at several frequencies. Generally, G{uj, Hq) 
depends on the frequency. Moreover, the real part of 



G {lu,Hq), changes its sign with increasing the dc field. 
The frequency dependence of G {oj, Hq) of the Nb crystal 
for several magnetic fields near the loss peak, presented 
in Fig. 1101 is shown in detail in Fig. Illb . For some mag- 
netic fields G {uj,Ho) changes its sign as the frequency 
increases. The zero in the real part of the response func- 
tion can be considered as the manifestation of the ex- 
istence of unknown collective mode in the system. In 
Fig. Illb we show the frequency of this hypothetical sur- 
face mode as a function of the dc field for Nb at 8.5K. 
The dynamics of K cannot be described by the first order 
linear differential equation: 



dK/dt = -v{K --ih{t)) (8) 

O.J(Hq.Ko) 

with some relaxation constant v and 7 = — aj{Hl°Ko) > 0- 

Equation ((HJ gives the real part G {uj,Ho) = "fu'^l{v^ + 
cj^), which does not change its sign with frequency. One 
could expect that the second order differential equation 
will give adequate description of the observed response. 
But our experimental data show, that G (lu^Hq) cannot 
be obtained from a differential equation of comparatively 
low order. 

The observed logarithmic frequency dependence of x j 
Fig. for some dc fields, resembles the response of a 
spin-glass system. But in spin-glass materials x^ is a 
slow function of the frequency as compared to x ■ Our 
data show that the x and x values, both depend on the 
frequency. Probably, similar to the spin-glass systems, 
we have here a lot of clusters which are governed by the 
second order differential equation, thus the observed re- 
sponse is the average over all clusters. The second order 
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FIG. 10: (Color online)The response function G{uj,Ho) of 
Nb and ZrBi2 samples (upper and lower panels, respectively) 
versus reduced magnetic field, Ho/Hc2- 
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FIG. 11; (Color online) (a) Frequency dependence of G at 
different magnetic fields near the absorbtion maximum for Nb 
crystal at T = 8.5 K. 

(b) Dispersion relation of the surface mode, u!o{Ho) — ^(y j P 
for Nb sample at T = 8.5 K 



differential equation 

dKldt = -v{K - fhit)) - f3d^K/dt^ (9) 

with /3 > has a response that changes its sign at 
u!q — v/p. Assuming that v/(3 for all clusters increases 
with the dc field, we find, that the dc field value for which 
G {uj,Hq) = must increase with the frequency, as ob- 
served in Fig. ^1 The \/v j fi quantity in this case will 
be the frequency of the surface collective mode. The dis- 
persion relation of this mode, loo{Hq), for Nb is shown in 
Fig.HIfc- 

Further insight into the low-frequency response can be 
obtained from the Kramers-Kronig relation [23| 



(10) 



where Xoo = Xs{oo). If we chose << ^ << 6m then: 

^ ,2( ria ix"ii)di roc Cx" (O^C N (11) 



The left side of this equation can be extracted from the 
available experimental data. By presenting it as a linear 
function of l/w^ one can obtain the averaged imaginary 
part at low frequency, defined as: 



?0 



(12) 



In Fig. 1121 we have shown the averaged imaginary part of 
G^(6o) (G" oc x'l), obtained for Nb at T = 8.5 K, using 
60 = 27r20 s-i,6,„ = 27rl065 s~\ G^(Co) in Fig. 12, 
is considerably larger than any value shown in Fig. 1101 



8 




1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 



H„/H , 

c2 



FIG. 12: (Color online) G" {lu, Ho) over the frequency interval 
< ujjl-K < 20 Hz as a function of reduced magnetic field, 
Ho/H^2, for Nb at 8.5 K and ZrBi2 at 4.5 K. 



So the Kramers-Kronig relations predict the existence of 
large loss peak at low frequencies. This prediction needs 
a further experimental evidence. 

The observed ac response of SSS has a very com- 
plex character. Partially, it can be the result of non- 
homogeneity of the samples. The K parameter, as it was 
introduced by Eq. corresponds to the wave function 
of the whole sample. In real samples, due to inhomogene- 
ity, K is not the integral constant of the GL equations 
and the exact wave function is the superposition of states 
with different K. Those states can relax with different 
relaxation times as in spin-glass systems ^3 j there- 
fore we observe such complicated responses. In addition, 
the dynamics of the surface state is not described by a 



first order simple relaxation equation. As a result, we 
obtain the logarithmic frequency dependence of the real 
part of the susceptibility, as in spin-glass systems, but the 
imaginary part shows a maximum at some frequency. 



VI. CONCLUSIONS 

In this paper we have presented investigation of the 
linear ac susceptibility of Nb and ZrBi2 single crystals 
in the surface superconducting state. Losses in this state 
have a linear origin, and the critical state model for the 
surface current does not apply here. Similar to spin- 
glass systems (where finite losses at considerably low fre- 
quencies exist), the real part of susceptibility exhibits a 
logarithmic frequency dependence. But the out-of-phase 
component has a frequency dispersion. This dispersion 
in SSS differs from that of the spin-glass systems. We 
assume that the sample surface presents a lot of super- 
conducting clusters, which are governed by second order 
differential equation, and the observed response is an av- 
erage over these clusters. The Kramers-Kronig analysis 
of experimental data reveals huge absorption peak at low 
frequencies. The response of K to the ac magnetic field 
defined by Eq. G{u!,Ho), has been measured. 
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